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We study the scattering properties of Sine Gordon kinks on obstructions in the form of finite size 
potential 'wells'. We model this by making the coefficient of the cos{(f) — 1 term in the Lagrangian 
position dependent. We show that when the kinks find themselves in the well they radiate and 
then interact with this radiation. As a result of this energy loss the kinks become trapped for 
small velocities while at higher velocities they are transmitted with a loss of energy. However, the 
' interaction with the radiation can produce 'unexpected' reflections by the well. 

■ We present two simple models which capture the gross features of this behaviour. Both involve 

, standing waves either at the edges of the well or in the well itself. 
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I. INTRODUCTION 



, Recently we have performed a detailed study of scattering properties of (2+1) dimensional topological solitons on 
, potential wells. This work was based on the 'baby' Skyrme model, ie we used the Lagrangian density which consisted 



of three terms: the pure S sigma model, the Skyrme and the potential terms: 



o . ^ ^ 

(-H ' where 

Oh! F(</.)=m(1 (2) 

D , ^ 

' The vector (j) was restricted to lie on a unit sphere hence we put (j) ■ (f) ~ 1 . 
^ ' To generate the potential well the coefficient of the potential term /j, was made x dependent. For x outside the well 
' it had one value, say, ^out and inside the well, ie for a < x < 6 its value was reduced to This choice of /i did not 
■ affect the vacuum (taken as 03 — +1); the skyrmion was given by a field configuration which varied from (p^ — +1 
far away from the position of the skyrmion to 03 = —1 at its position). Initially, the skyrmions were placed far away 
from the 'well' (so that all the variation of 03 from +1 took place for x well away from the 'well', ie for x < a). 

The 'skyrmions' were then sent towards the well and their properties were studied. The obtained results have shown 
that when the solitons fall into the 'well' they get deformed and this deformation may excite the vibrational modes 
of the skyrmion and may lead to the skyrmions radiating away their access of energy. In consequence, the skyrmions 
can get 'trapped' in the well or emerge from it with a reduced velocity. In Q we presented a simple four mass model 
which apes the vibrational modes of the skyrmion and we showed that many of the observed scattering properties 
of skyrmions can be reproduced in this model - suggesting that their origin resides in the excitation of the lowest 
vibrational modes of the skyrmion. 

Given this, it is important to check what happens in models in which solitons have fewer vibrational modes and so 
we have decided to look at the (1+1) dimensional Sine-Gordon model and see what happens when its kinks scatter 
on the potential wells. 

In the next section we discuss our results obtained for the Sine-Gordon model. In this model we include a finite 
size, finite depth, potential well which is introduced by appropriately modifying the coefficient of the nonlinear term 
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in the Lagrangian. The results are quahtatively similar to what we have seen in the two-dimentional model and are 
not very different from the results obtained some time ago by Fei et al fsf , in a work which involved the scattering of 
Sine-Gordon kinks on a one-point impurity^. 

Recently, Goodman and collaborators [4] have explained these old results jSj in a 'two bounce' resonance model. 
Their explanation is based on the interaction of the kink with the oscillation of the vacuum (around the impurity). 
Thus their model involves two degrees of freedom - the position of the kink xq and the amplitude of the vacuum 
oscillation (at the impurity point) a. The model of Goodman et al has reproduced all the features of the results of the 
original simulations reported in |3|. Hence in the following section we introduce a similar model for our case which 
now involves a finite well of width 2p and depth 1 — A. To do this we make an ansatz for an approximate field which 
describes the system. It involves a sine-Gordon kink which is allowed to alter its slope and we add to it two amplitudes 
of oscillation of the vacuum (at each end of the 'well'). In section 3 we derive the Lagrangian for such an effective 
model from the original Lagrangian. As the model is somewhat crude we make some drastic approximations in our 
derivation of the Lagrangian but still find that the model reproduces the main features of the scattering reasonably 
well. Hence we believe the ideas of Goodman and collaborators to be correct and be more general in nature - thus 
showing that due to the interaction of the soliton with the radiation in the well its behaviour can be quite complicated 
and can result in the reflection of the soliton by the well; ie a process which is purely classical in nature but could be 
confused with a quantum behaviour. To confirm this further we introduce a further model (with a couple of radiation 
standing waves in the well) and again reproduce the main features of the full scattering process. 

The last section presents some concluding remarks. 

II. SINE GORDON MODEL AND ITS KINKS. 

We take the Lagrangian of the (l-fl) dimensional Sine-Gordon model in the form 

£ = • - sin((^)^ (3) 

where, for the kink, the basic field (p goes from at a; = — oo to tt at = oo. In the static case its explicit form is 

tp{x) = 2tan"-^(exp(6'(a; - .To)), (4) 

where xq is the kink's position and 9 is its slope. For ([4| to be a solution of the equations of motion which follow 
from ([3]) we need to set 6* = A. 
To have a 'well' we set 

{1 for \x\ > p 
(5) 
Ao for < p. 

Clearly, Aq < 1 describes a well while Ao > 1 describes a barrier. As the two-dimensional studies gave more interesting 
results for the wells in this paper we restrict our attention to Aq < 1. 

We have performed many numerical simulations, varying Aq and p (the width of the well). 



A. Numerical Simulations 



We have performed most of our simulations using a 10001 point lattice with the lattice spacing being 0.05. Hence 
our lattice extended from -50 to +50. The kink was initially placed at xq = —40. Its size was determined by 6* = 1 
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which means that its field was essentially ip ^ ior x < —45 and ~ tt for x > —35. Thus, there were no problems 
with any boundary effects (wc have verified this by altering the lattice size and xq). 

We have performed three sets of simulations; one involving a very narrow 'well' {p = 0.5) and two involving a larger 
well {p = 5) (one shallow - Ao = 0.8 and one rather deep - Ao = 0.2). In each case we sent the kink (originally at 
xo = —40) towards the well varying its initial velocity. To do this we took the expression for a kink moving with 
velocity v ie 

'fi{x,t) = 2 ta.n~ ^{exp{'-f{x — Xo — vt j), 7 = — ^ (6) 

Vl — 

(obtained by Lorentz boosting (4) and setting ^ = 1) and then used it to calculate the initial conditions (<^(a;,0) and 
at /• 

All our simulations were performed using a 4th order Runge-Kutta method for simulating the time evolution. 
The time step of our simulations was taken to be 0.0001. We used fixed boundary conditions and later we used 
also absorbing boundary conditions at the edges of the lattice. This was generated by successively decreasing the 
magnitude of ^ for the last 50 points at both ends of the lattice. 

B. Deep Well 

In the deep well case we took Aq = 0.2 and p ~ 5. We have found that when the kink was in the well it radiated 
and so, when it finally emerged from the well its velocity was lower than the initial velocity. This was due to the fact 
that the well distorted the soliton which then began vibrating (ie its slope started oscillating). These vibrations were 
then gradually converted into radiation with the slope settling at its original value. The curve of the outgoing velocity 
as a function of the incoming one is shown in fig. 1. As is clear from the plot, the kink whose initial velocity is less 
than Vcr ^ 0.527 gets trapped in the well. The curve is very smooth and, as expected, we note that for incoming 
velocities larger than Vcr the outgoing velocity is always larger than the incoming one demonstrating the loss of the 
kinetic energy of the kink through vibration resulting in radiation. 

C. Shallow Well 

For a shallow well we took Aq = 0.8 and still p = 5. This time the critical velocity is much smaller (as the well 
perturbs the kink much less). The curve of the outgoing velocity as a function of the incoming one is shown in fig. 
2a. 

Looking at the plot we note that there is some irregularity close to the threshold. Blowing it up (see fig. 2b) 

wc note that just below the critical velocity we also have some negative velocities (ie the whole process looks as if 

the kink was reflected by the well!). Thus in addition to trajectories like those in fig. 3a and fig. 3b we also have 

trajectories like those of fig. 3c 

Of course - the reflected trajectory is somewhat unexpected; this is what one would expect in a quantum system 

but here we have a completely classical system and we have a reflection by the well. Clearly, this reflection must be 

somewhat related to the interaction of the deformed kink with the radiation that is present in the well. To look at 

this in more detail we have looked also at a very narrow well. 
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FIG. 1: Velocity of the kink leaving the well as a function of its velocity as it approaches the well, p = 5, Ao = 0.2. 

D. Very Narrow Well 

This time we took p — 0.5 and Aq = 0.2. The obtained curve of the velocities is shown in fig. 4a. Again we see an 
interesting behaviour close to the threshold - and blowing it up we get fig. 4b. 

Are these results just numerical artifacts? To answer this question we have carried out several tests. First we 
changed the lattice size (increased the number of points, changed the lattice spacing) and changed the time step in 
the Runge Kutta scheme. However, the observed pattern of reflections was reproduced in all simulations. We have 
also varied the discretisation scheme and, among others, considered also the 'topological discretisation' of Speight 

n 

and Ward [5;]. All simulations using these discretisations exibited similar patterns, with the values of the velocities 
esentially unchanged (the values changed in the fourth decimal points) . Hence we believe the effect to be genuine and 
so we are left with having to explain its origin. 

Performing a literature search we have found the above mentioned paper by Fei et al . In that paper the authors 
studied the scattering of Sine-Gordon kinks on a one lattice point impurity (at one lattice site the potential Xtp'^ was 
changed to a different value - X'lp'^). Fig. 11 of that paper looks amazingly similar to our fig. 4b. Of course, our 
potential corresponds to a superposition of defects of Fei et al so not surprisingly the pattern is more evident in a 
system involving a smaller well. However, the fact that the effect persists and is seen for larger wells suggests that 
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FIG. 2: Velocity of the kink leaving the well as a function of its velocity as it approaches the well {p = 5, Ao = 0.8). a) Full 
plot, b) A close up for small velocities. 



the phenomenon is more fundamental in nature. 

The phenomenon observed by Fei et al was recently explained by Goodwin et al in terms of a two mode model 
invoving the interaction of the kink with the radiation at the impurity point. Clearly in our case we have two special 
points (at the two 'edges' of the well where the potential changes). Hence it may make sense to develop a model based 
on osciallations at these two points. This is what wc do in the next section. 



III. OUR EFFECTIVE MODEL 

Following Goodman et al Q we consider the following ansatz 

(pix,t) = 2 tan"^ exp{9{t){x ~ xo{t))) + a{t)g{x) + b{t)h{x), (7) 
where we take g{x) and h(x) in the form: 

g{x) = exp (~^-^d] , h(x) = exp (~^-^f] , (8) 





and 

D= - ^ , F ^ l"' - ^ . (9) 

\B,x<p \A,x< 

Thus we have allowed the kink to change its position xo{t) and its slope 0{t), and a{t) and b{t) represent the excitations 
of the vacuum. 

To obtain an effective model we put this expression for (p{x,t) into the Lagrangian density ([3]) and attempted to 
integrate it over x. However, as this leads to rather complicated expressions which we had not succeeded to calculate 
analytically we resorted to an expansion (in the pcturbation due to the well). Hence we put ip — tpo + Sip, where 

ipo = 2 tan^i exp i0{t){x - xo{t))) (10) 
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FIG. 3: Typical trajectories (p = 5, A = 0.8). a) Below 'threshold' {v = 0.071535), b) Transmitted {v = 0.077983), c) Reflected 
{v = 0.071514). 



and 



S<fi = a{t)g{x) + b{t)h{x). 



(11) 



Next we expanded: 



s\v? [ifQ + 5^p) ~ sin((^o) + sin(2(po) '^V' + cos{2ip) {5ipY ... 



(12) 



Then as 



6{x-xo) - 0X0 e^(^-^«) + dg{x) + bh{x) 



(13) 
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FIG. 4: Velocity of the kink leaving the well as a function of its velocity as it approaches the well {p = 0.5, Ao = 0.2). a) Full 
plot, b) A close up for small velocities. 



we found that 



The 6ip terms give us by themselves 

1 



J — oo 



dt 



ah 



12 6*3 



-Bp 



+ 9x1 



-pB 



-pA 



A+B A-B^ ' 

Finally, we calculated the 'crossed terms' and we got (for the terms involving xq and d and 6) 



-26xq a 



K 



1 



1 



~ 2 2 



where K = exp{9{p — xq) and a similar expression for —26xob, except that this K was replaced by 
K = exp(^(— p — a;o)). For the terms involving 6 and a we got 



20 d 



K 



1 + 



1 



1 



(14) 



(15) 



(16) 



(17) 



The expression involving Ob was again the same with K replaced by K. These expressions are not exact; we obtained 
them by making several approximations of the type 



K 



,-(0 + ^ 



K 



1 



+ K'^e-^^^ 1 + (6" + f ) ' 



(18) 



which should be reasonably reliable given the exponential form of the dependence of the integrands on z etc. 

Next we calculated the contribution due to (§f Performing the integrations, and making similar approximations 
as before we got 

(19) 



8 4 



4(A + B)e-^- + ^^±^(e-^-e-^) 
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plus terms linear in a and h. They are given by 



^-0 f +. 



(20) 



Finally we had to add the A sin (ip) terms. The contribution due to the well is given by 



Aq / sm^{<po)dx = Aq y 



P g2e(x-xo) 



1^1 -I- g2e(x-xo^2 



dx 



(21) 



_ Ag sinh(2gp) 

~ T (cosh(26ia;o) + cosh(26'p) ' 

The contribution due to dx Ag sin^((^) can be calculated in a similar way and we get 



P. = 4 



(22) 



from the (fo term, 



K{l-K^) K{l-K') 



(23) 



from the terms linear in 6(p and 



P5 



A? 



A + B 



-Bp 



A-B 



-pB _ --pA> 



(24) 



Adding all these terms together we get the Lagrangian for our effective model. This Lagrangian is given by 
i = ^ J + ^i;' + lia" + b^)Bo + abBi - 2exo{aDo{6) + bD^{6)) 



20{dD2{9) + bD3{0)) 



^^^iA + B)-^^ + ^^eiaDoiO)+bD,ie)) 



(25) 



A? 



ab . 



2eXlDe{9,xo) - -f - 2Xl{aDi{e) +bD^{6)) - Xf-B^ - Xf 



0) 



where 
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Bi 



-pB 



A + B A-B 



-pB _ ^-pA\ 



and where the 7 frmctions Di i = 0,..6 are given by: 



Do{0) = 



K 



1 + ^2 



1 



1 



' + 0.5^-61 



K 



l + K'^ 



r) + 



(61 + 0.55)2' (0.5A-6i)2 



+ P 



1 + 



1 



+ 



1 



6I + 0.5B 0.5A-6I 



D4i6) = Do{e) 



1 + 
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The functions Di{0) and D^(0) have the same form as Dq(9) and Di{9), respectively, after the replacement K ~* K 
and 



DM = 



K 



Finally Dq{9,xo) is given by 



1 



O.dB)' 



r) + 



1 



1-K 



2 r 



1 



iO.5A-0)^ 

sinh(26'p) 



1 



1 



9 + 0.5B 0.5A-9 



9 cosh(6'(a;o + p)) cosh{9{p — Xq)) ' 
The derived Lagrangian involves 4 variables XQ,9^a and b but its equations are rather complicated. Hence we 
have solved them numerically. In the next section we discuss our solutions of the equations which follow from the 
Lagrangian ([25]). 



IV. RESULTS IN OUR EFFECTIVE MODEL 



We have looked at the equations for 9, xq and a and b, which follow from the Lagrangian ([^5]) and solved them 
numerically. As our initial conditions we took xq — —40, 9—1 and set, initially, a = b = 0. Of course, we also put 
a = b = 9 = and studied the behiour of our system as a function of xq. We also varied a little the parameters A and 
B, which appear in the description of the effects due to the well. For most of our work we used the values around 
0.5 (and used the fact that we expect A^ ^ e + B^. To have the well similar to the one we used in the full simulation 
(shallow well) we put e ~ 0.16. This is due to the fact that the linearised equations (ze for the waves) differ, inside 
and outside the well, by a term proportional to (Aq — l)u. In our case this translates to A^ — B^ ^ e and so for a 
shallow well we have e ^ 1 — Aq = 0.16. We simulated the time evolution using the 4th order Runge Kutta method 
and have found that the well distorts the kink quite strongly but, as expected, it can trap the kink like in the original 
Sine-Gordon model. However, as the effective system does not absorb energy after a few bounces the kink can escape 
(forwards or backwards) . Of course in the real system there are many degrees of freedom of radiation which can get 
generated in the hole and such bounces are very rare. So to model these 'extra' modes of radiation which take the 
energy out of the modes we are describing we introduced an absorption of the oscillations of 9{t), ie. we added a 
term proportional to 9 in the equation for 9. This has, indeed, reduced the oscillations in 9 and made the model more 
realistic. 

In fig. 5 we present the curve of Vout as a function of Vin obtained in our model, ie based on the Lagrangian (j25p . 
We note that although our effective model is quite crude it reproduces fig. 3b rather well. Hence we believe that the 
mechanism of Goodman et al P| explains the behaviour of the kinks in our case too. Of course, we could check this 
in more detail by making fewer approximations and, perhaps, performing the evaluation of all the difficult integrals 
numerically, but we are not sure that the extra effort required would be justified. 



V. ANOTHER EFFECTIVE MODEL 



Given that our effective model reproduces the results of our simulations rather well we have decided to check the 
generality of this observation - is this related to the existence of standing waves in tbe neighbourhood of the well or 
is our model somewhat unique. Hence we considered another model this time involving a standing wave in the well. 
The idea here is to explore further whether the relections observed in the full model can be related to the interaction 
of the soliton with the waves in the well. Hence this time we have taken for our field configuration ip — ipo + Sip, 
where, as before, 

ipo = 2 tan^i exp i9{t){x - xo{t))) (26) 
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FIG. 5: Outgoing velocities as a function of incoming ones - in our effective model. 



and 



This time for g{x) we have taken 



6^ = a{t)g{x). (27) 



5(2^) = cos(— ) - -cos(— ) (28) 

inside the well {ie for |a;| < p and zero otherwise). Our g{x), and its derivative, vanish at x = ztp and the idea is that 
as the soliton enters the well its slope 9{t) changes. This excites the modes described by g{x) and so a{t) becomes 
nonzero. This puts some energy into these modes which then interact with the soliton. Due to this interaction this 
energy can, every now and then, be given back to the soliton resulting in its reflection or transmission. 

Like in the previous model, we have put the expression for (p into the original Lagrangian, integrated over x obtaining 
an effective model involving xq, 6 and a. Then we performed some simulations starting with a soliton originally far 
from the well moving with some velocity towards it. Like in the previous model we have absorbed energy through a 
term propertional to ^ in the equation for 6. Also, like in the previous case wc have found some reflections below 
the threshold for the scattering through the well. Their details depend a little on the strength of the absorption; at 
no absorption we have reflections, and transmissions forward, at much lower velocities, the higher absorptions kill 
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FIG. 6: Outgoing velocities as a function of incoming ones - in our second effective model. 



the reflections and transmissions. They also raise the threshold for transmission. Of course, it is difBcult to estimate 
reliably the degree of absorption as in the full model there are many modes of oscillations in the well. So in the end 
we have used very low absorption; then the threshold for transmission was close enough to what is seen in the full 
simulation but we had more reflections below threshold. To improve the model we should have used more modes in 
the well but this would have resulted in a more complicated model with more degrees of freedom. As wc have only 
wanted to test our basic idea (that the reflections are due to the interaction of the soliton with radiation modes), 
in this paper, we decided to look at a model which is as simple as possible ie with a very few degrees of freedom. 
Our choice of two low lying modes was dictated by simplicity and the basic belief that the lower modes would get 
excited first and so are more important. In fig. 6 we present the curve of Vout as a function of Vin obtained in this 
(second) model (with a reasonble absorption) . Once again, we see that the model reproduces well the pattern of the 
simulations of the full Sine-Gordon model (fig. 3b). 



VI. CONCLUSIONS 



We have looked at a system involving a Sine-Gordon kink scattering on a 'well'-like potential. 

We have found that, like in (2+1) dimensions, when the kink was sent towards the well it gained some energy as it 
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FIG. 7: Outgoing velocities as a function of incoming ones - in our second effective model without relativistic corrections. 



entered the well. Some of this energy was converted into kinetic energy of the kink, some was radiated away. So when 
the kink tried to 'get out' of the well it had less kinetic energy than at its entry and, when this energy was too low it 
remained trapped in the well. However, as the soliton moved in the 'well' it interacted with the radiation in the well 
and at some specific values of the initial velocity this interaction resulted in the kink being ejected backwards from 
the well (with much reduced velocity). Thus, seen from outside, the well acted as if it reflected the kink, something 
which is seen in quantum systems but which is less well known in classical systems. 

We have performed many numerical simulations to make sure that the observed behaviour is not an unexpected 
artifact of our numerical procedures and the pattern survived all applied tests. Hence we believe the effect to be 
genuine. 

We have noted that a similar behaviour was observed many years ago by Fei et al who studied the scattering of 
kinks on a one-point impurity. This behaviour was recently explained by Goodman et al Q as a two-bounce resonance 
between the kink and the oscillation of the defect. This has made us to consider two models of a similar nature. Both 
models are very simple, clearly too simple, but we wondered whether they would qualitatively reproduce the observed 
effects. Both involve kinks interacting with radiation and we generate this by taking an ansatz involving a kink (which 
can vary its position and slope) plus a couple of radiation modes. This ansatz is then put into the full equation which 
are then integrated out resulting in a model involving a few variables; namely, the position and the slope of the kink 
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FIG. 8: Outgoing velocities as a function of incoming ones - in the full simulation without relativistic corrections 



and the coefficients of the oscillation modes of the vacuum (modeling radiation) . 

The first model involved taking standing waves that are located at the eges of the well; the other one involved just 
one wave in the well. We have found that both models reproduced the main features of the observed behaviour quite 
well suggesting that the mechanism of Goodman et al is more general in nature and that, in general, the observed 
phenomenon of reflection of the solitons on the well is related to their interaction with the waves in the well. Of course, 
both models are too simplistic; to describe properly the full process we have to understand better which modes of 
radiation are important and why. This involves more work and is planned for the future. However, the work done so 
far suggests that we are on the right track and gives us encouragement for the further study. This is confirmed further 
by what we saw in a two-dimensional model [2| in which the solitons have genuine vibrational as well as radiation 
modes. 

Incidentally, in our calculation we have used the correct initial condition ie with the correct relativistic factors 
(--y = --^=i=). Had we ignored them and used their norelativistic form {ie not modified 9 in (10) and (26)) we would 
have obtained instead of Fig. 6 the dependence which is shown in Fig. 7. We note an interesting oscillatory behaviour. 
Of course, this oscillation is overemphasised by the use of too few radiation modes but we wondered whether it would 
be seen in the full model too {ie whether the addition of further modes would wash them out). Hence we have redone 
the full simulations also without the relativistic factors. Our nonrelativisitc curve is shown in fig. 8. We note the 
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extra oscillations. Their origin lies clearly in the fact that the absence of 7 factors induces initial distortions resulting 
in the change of 6. This affects the phase of the soliton and so alters its interaction with the wave in the well resulting 
in the curve shown in fig. 8. It is interesting to note that our effective models also reproduce these oscillations, thus 
giving further support for the validity of our claim that the interaction with the well proceeds through the generation 
of standing waves and their intereference with the solitonic fields. Of course the agreement between the results of 
the full simulations and of the efii'ective models is only qualitative in nature as in our effective models we used only 
some standing waves which were chosen somewhat ad hoc. To get a quantitative agreement we have to determine the 
relative importance of different waves - this problem is CTirrcntly under consideration. 

Several real physical systems arc described by the sine Gordon equation, especially in solid state physics, and it 
would be interesting to see what the physical implications of a position dependant potential, like the one used in this 
paper, would be. 
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